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ABSTRACT 

Using a large sample of 38,478 star-forming galaxies selected from the Second Data Release of 
the Sloan Digital Sky Survey database (SDSS-DR2), we derive analytical calibrations for oxygen 
abundances from several metallicity-sensitive emission- line ratios: [N nJ/Ha, [O IIl]/[N n], [N n]/[0 n], 
[N n]/[S n], [S n]/Ha, and [O iii]/H/3. This consistent set of strong- line oxygen abundance calibrations 
will be useful for future abundance studies. Among these calibrations, [N n]/[0 n] is the best for 
metal-rich galaxies due to its independence on ionization parameter and low scatter. Dust extinction 
must be considered properly at first. These calibrations are more suitable for metal-rich galaxies 
(8.4<12+log(0/H)<9.3), and for the nuclear regions of galaxies. The observed relations are consistent 
with those expected from the photoionization models of Kewley & Dopita (2002). However, most of 
the observational data spread in a range of ionization parameter q from 1 x 10 7 to 8 x 10 7 cm 
s _1 , corresponding to \ogU~~ 3.5 to —2.5, narrower than that suggested by the models. We also 
estimate the (N/O) abundance ratios of this large sample of galaxies, and these are consistent with 
the combination of a "primary" and a dominant "secondary" components of nitrogen. 
Subject headings: galaxies: abundances - galaxies: evolution - galaxies: ISM - galaxies: spiral- 
galaxies: starburst 



1. INTRODUCTION 

The chemical properties of stars and gas within a 
galaxy are a fossil record chronicling its history of star 
formation. Optical emission lines from H n regions have 
long been the primary means of gas-phase chemical di- 
agnostics in galaxies (Aller 1942; Searle 1971; reviews by 
Peimbert 1975; Pagel 1986; Shields 1990; Aller 1990; Pe- 
imbert, Rayo, & Torres-Peimbert 1975; Osterbrock 1989; 
de Robertis 1987). Osterbrock (1989) thoroughly dis- 
cusses the standard techniques for measuring the chemi- 
cal properties of ionized gas, which require measurements 
of H and He recombination lines along with collisionally 
excited lines from one or more ionization states of heavy 
element species. 

Oxygen is the most commonly used metallicity indi- 
cator in the interstellar medium (ISM) by virtue of its 
relatively high abundance and strong emission lines in 
the optical region of the spectrum (e.g., [O n] A 3727 
and [O iii]AA4959, 5007). Ideally, the oxygen abundance 
is measured directly from ionic abundances obtained 
through a determination of the electron temperature (T e ) 
of the H II region. This is known as the "T e — method" 
(Pagel et al. 1992; Skillman & Kennicutt 1993). The ra- 
tio of a high-excitation auroral line such as [O iii]A4363 
to the lower excitation [O iii]AA4959,5007 lines provides 
a direct measurement of the electron temperature in the 
physical medium where ++ is the dominant ionization 
state. However, [O iii]A4363 is only prominent in metal- 
poor H II regions. In most cases, where [O iii]A4363 is 
too weak to be detected, metallicities of H n regions have 
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to be estimated from "strong-line" ratios. A commonly 
used one is R23: 



R23 = 



/([OIII]AA4959, 5007) + /([OII]AA3726, 3729) 

JW) 



(1) 

Many calibrations of converting R23 to 12+log(0/H) 
abundance have been published, including Pagel et 
al. (1979), Pagel, Edmunds, & Smith (1980), Ed- 
munds & Pagel (1984), McCall, Rybski, & Shields 
(1985), Dopita & Evans (1986), Torres-Peimbert, Peim- 
bert, & Fierro (1989), Skillman, Kennicutt, & Hodge 
(1989), McGaugh (1991), Zaritsky, Kennicutt, & Huchra 
(1994), Kobulnicky et al. (1999, K99), Pilyugin 
(2000, 2001a,b), Tremonti et al. (2004, T04), Salzer 
et al. (2005) and Kewley & Dopita (2002) (here- 
after KD02) etc. Some calibrations have attempted to 
take the ionization parameters into account (e.g., Mc- 
Gaugh 1991; KD02 etc.), but others have not. How- 
ever, one problem is that the relationship between R23 
and 12+log(0/H) is double valued, with the transi- 
tion between the upper metal-rich branch and the lower 
metal-poor branch occurring near 12+log(0/H)^8.4 
(logi?23^0.8). Moreover, in many cases R23 is unavail- 
able due to the limit of wavelength coverage or the 
poor quality around one or some of the related lines. 
One can then use other various metallicity-sensitive 
strong-line ratios to estimate metallicities of galax- 
ies, including [N n]A6583/Ha, [O m]A5007/[N n]A6583, 
[N n]A6583/[0 n]A3727, [N n]A6583/[S n]AA6717,6731, 
[S n]AA6717,6731/Ha, and [O ni]AA4959,5007/H/3. Even 
when R23 is available, these ratios are also useful for 
breaking the R23 upper/lower branch degeneracy. We 
should notice that some of these strong-line ratios are 
also sensitive to interstellar reddening (e.g. N 11] /[O 11]) 
or depend on the ionization parameter. 
Some researchers have used T P -based abundance data 
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to calibrate for some of these strong-line ratios, for exam- 
ple, Denicolo et al. (2002) obtained a linear calibration 
of [N n]/Ha to 12+log(0/H) from a combined sample 
of 108 metal-poor galaxies with [O in] 4363 detected and 
128 metal-rich galaxies; Pettini & Pagel (2004) studied 
the calibration relations between [N n]/Ha, ([O m]/H/3 
)/([N il] /Ha) and 12+log(0/H) from 137 galaxies with 
T e -derived O /H abundances (including six of them from 
photoionization models). Perez-Montero & Diaz (2005) 
also studied some calibrations from the collected 367 
emission-line objects on the basis of the T e -based oxygen 
abundances, but they did not derive analytic calibrat- 
ing relations. However, these observations are limited to 
modest samples, and only apply to some of the strong- 
line ratios. The much larger consistent data set (sev- 
eral tens of thousands) from the recently released SDSS 
database could provide much more information about 
many of the strong-line ratio calibrations for metallici- 
ties of galaxies. 

KD02 used a combination of stellar population syn- 
thesis and photoionization models to calibrate several 
strong-line abundance diagnostics. In the KD02 models, 
strong line ratios are a function of both metallicity and 
ionization parameter. The ionization parameter, q, is a 
measure of the number of ionizing photons per atom at 
the inner boundary of the nebula. To determine both 
metallicity and ionization parameter, KD02 suggested 
a somewhat complicated iterative approach involving a 
number of line ratios. We advocate a simpler approach in 
this study. In Sect. 4 we compare the observed line ratios 
of a large sample of star-forming galaxies selected from 
the SDSS-DR2 to the predictions of the KD02 models 
and show that galaxies span a smaller range of q than 
used by KD02. We therefore ignore any explicit depen- 
dence of the line ratios on q, and use observational data 
to derive strong-line abundance calibrations which reflect 
the mean ionization parameter of the data at a given 
metallicity. 

The SDSS-DR2 database is very well suited for this 
kind of study. There are a large number of galaxies at 
any given value of oxygen abundance between about 8.4 
and 9.3, which makes it ideal for measuring correlations 
between metallicity and various strong-line ratios and 
for quantifying the scatter. We make use of the emis- 
sion line data and metallicities described in Brinchmann 
et al. (2004) and Tremonti et al. (2004) which have 
been made publicly available on the web 4 . The oxygen 
abundances presented in Tremonti et al. (2004) were 
derived using Bayesian techniques to compare multiple 
strong emission lines to a grid of photoionization mod- 
els (Chariot et al. 2006). We use the R23 metallicity 
calibration of Tremonti et al. to derive a fiducial oxy- 
gen abundance. We then examine correlations between 
the oxygen abundances and the observed several strong 
emission- line ratios: [N 11] /Ha, [O m]/[N 11], [N 11] /[O 11], 
[N 11 ]/[S 11], [S 11] /Ha, and [O m]/H/3, then derive 
new strong-line metallicity calibrations where appropri- 
ate. Our approach differs from the usual way in which 
abundance indicators are calibrated: 1) using samples of 
galaxies with direct (T e -based) abundances (e.g. Pagel et 
al. 1979; Pettini & Pagel 2004); 2) using photoionization 
model grids (e.g. McGaugh 1991; Chariot & Longhetti 

4 http://www.mpa-garching.mpg.de/SDSS/ 



2001; KD02). Our calibrations are valid for high metal- 
licity (12+log(0/H) > 8.4) galaxies, where the T e -based 
abundances arc only available for a handful of objects 
(Bresolin et al. 2004, 2005; Garnett et al. 2004a,b)! We 
note that our calibrations may be subject to the same 
systematic uncertainty as the models. For comparison, 
we use the R23 calibration of McGaugh (1991) (as Kob- 
ulnicky et al. 1999) to derive the fiducial metallicities as 
well. 

We also estimate the log(N/0) abundance ratios of 
these sample galaxies by using the algorithm given by 
Thurston et al. (1996), including the estimated temper- 
ature in the [N 11] emission regions. This will be the first 
such large sample for the observed log(N/0) abundances, 
and will be useful for understanding the "primary" or 
"secondary" origin of nitrogen. Specially, the strength 
of the SDSS dataset is that there are a large number of 
galaxies at any given value of oxygen abundance between 
about 8.4 and 9.3. This makes it ideal for quantifying the 
scatter in various trends. 

This paper is organized as follows: the sample selec- 
tion and R23 -derived oxygen abundances are described 
in Sect. 2; Sect. 3 shows the derived calibrations of the 
linear and/or third-order polynomial fits from the ob- 
servational relations of 12+log(0/H) versus the various 
strong emission-line ratios; these observational relations 
are compared with the photoionization models of KD02 
in Sect. 4; The log(N/0) abundance ratios of the sample 
galaxies are given in Sect. 5; and we conclude in Sect. 6. 

2. SAMPLE SELECTION, AND R23-DERIVED OXYGEN 
ABUNDANCES 

The data analyzed in this study are drawn from the 
SDSS-DR2 (Abazajian et al. 2004). These galaxies are 
part of the SDSS "main" galaxy sample used for large- 
scale structure studies (Strauss et al. 2002). We select 
the sample galaxies having Petrosian r magnitudes in the 
range 14. 5< r < 17.77 mag after correction for foreground 
Galactic extinction using the reddening maps of Schlegcl 
et al. (1998), which yields 193,890 objects from the total 
of 261,054 objects. 

We only consider galaxies whose metallicities have 
been estimated by Tremonti et al. (2004), leaving 50,385 
galaxies. They selected star-forming galaxies using the 
criteria given by Kauffmann et al. (2003b) following 
the traditional line diagnostic diagram of [N Il]/Ha vs. 
[O iii]/H/3 which has been used by Baldwin, Phillips, Ter- 
levich (1981, BPT), Veilleux & Osterbrock (1987) and 
Kewley et al. (2001) etc. 

Tremonti et al. (2004) discussed the weak effect of 
the 3" aperture of the SDSS spectroscopy on estimated 
metallicities of the sample galaxies with redshifts 0.03 < 
z < 0.25. Kewley et al. (2005) further recommended 
that, to get reliable metallicities, redshifts z > 0.04 are 
required for the SDSS galaxies to ensure a covering frac- 
tion >20% of the galaxy light. Thus, we select the galax- 
ies with redshifts 0.04 < z < 0.25, which then leaves 
40,693 objects. 

We only consider the objects for which fluxes have been 
measured for [O 11], H/3, [O ill], Ha, [N 11] (also [S 11] 
in the study in Sects. 3.5, 3.6). It leaves 40,293 objects 
(39,919 if the [S n]6717,6731 are included to be consid- 
ered as well) . Also, the sample objects should have higher 
S/N ratios by requiring them to have lines of H/3, H a, 
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and [N n]6583 detected at greater than 5a. It leaves 
39,029 objects. 

The oxygen abundances of the sample galaxies were es- 
timated using the R23 method. The emission-line fluxes 
must first be corrected for dust extinction. The extinc- 
tion correction of the sample galaxies are derived using 
the Balmer line ratio Ha/H/3: assuming case B recombi- 
nation, with a density of 100 cm -3 and a temperature of 
10 4 K, and the predicted intrinsic ratio of Ha/H/3 is 2.86 
(Osterbrock 1989), with the relation of 

(h±) obs = ( ^) intr lo-(/(^)-/(^)). (2) 

Using the average interstellar extinction law given by Os- 
terbrock (1989), we have /(Ha) - /(H/?)=-0.37. Then, 
the extinction parameter Ay is calculated following 
Seaton (1979): A v = E(B - V)R = (mag). R= 
3.1 is the ratio of the total to the selective extinction 
at V. The emission-line fluxes of the sample galaxies 
have been corrected for this extinction. For the 311 data 
points with Ay<0, we assume their Ay=0 since their 
intrinsic Ha/H(3 may be lower than 2.86. 

Tremonti et al. (2004) used the approach outlined by 
Chariot et al. (2006), which consists of estimating metal- 
licities statistically, based on simultaneous fits of all the 
most prominent emission lines ([O 11], H/3, [O in], He 1, 
[O 1], Ha, [N 11], [S 11]) with a model designed for the 
interpretation of integrated galaxy spectra. Here, we use 
the corresponding emission-line fluxes to directly get dust 
extinction-corrected R23 values first, and then adopt the 
calibration given by Tremonti et al. (2004) (their eq.l: 
12 + log(0/H) = 9.185 - 0.313a; - 0.264x 2 - 0.321a; 3 , 
with x — logi?23) to estimate the 12+log(0/H) abun- 
dances for our sample galaxies, which we denote as 
12±log(0/H)^ 23 . This calibration is only appropriate 
to the metal-rich branch. The consistency of these two 
estimates is shown in Fig. Q] The solid line is the line 
of equality, and the two dashed lines are ±0.15 dex dis- 
crepancies from the line of equality. We notice that, for 
about 260 data points with 12 + log(0/H) < 8.5 esti- 
mated by Tremonti et al. (2004), our i?23-derived abun- 
dances are higher than those of Tremonti. All of these 
galaxies have -1.2<[N n]/Ha<-0.8 (see Fig. Ot), which 
implies that they are metal poor galaxies on the lower 
i?23 branch. Because the calibration of Tremonti et al. is 
only suitable for the upper branch, we select the galaxies 
having 8.4<12±log(O/H)<10 (39,006 objects). We also 
require the two metallicity estimates to be consistent to 
±0.15 dex, which leaves a final sample of 38,478 galaxies. 

If a histogram distribution of the 12±log(0/H)fl 23 is 
presented, it will show a range of 8.4<12 + log(0/H) 
<9.3 with a peak around 12±log(O/H)=9.0 and an in- 
creasing gradient-like distribution following the increas- 
ing abundance from 8.4 to 9.1. Among these selec- 
tion criteria, only the 0.04 lower limit of redshift, the 
8.4<12±log(O/H)<10 and the 0.15 dex discrepancy were 
not considered by Tremonti et al. (2004), others are the 
same. 

3. DERIVING CALIBRATIONS FOR OXYGEN 
ABUNDANCES FROM THE OBSERVATIONS OF THE 
SDSS GALAXIES 

We derive the analytic linear least squares and/or 
third-order polynomial fits from the observed relations of 



12±log(0/H)^ 23 versus "strong-line" ratios of the SDSS 
sample galaxies: [N 11] /Ha, [O m]/[N 11], [N n]/[0 11], 
[N 11] /[S II], [S 11] /Ha, and [O III] /H/3. To minimize the 
systematic uncertainty from calibrations of i?23 to O/H, 
we re-derive these calibrations by using two other sets of 
oxygen abundance estimates replacing 12±log(0/H)fl 23 : 
12±log(0/H)x99, the i?23-derived abundances by using 
the -R23-O/H calibration of Kobulnicky et al. (1999, from 
McGaugh 1991), and 12+log(O/H) T04 , the Bayesian 
abundances given by Tremonti et al. (2004). All the cal- 
ibration coefficients and the rms derivations have been 
given in Table 1. 12±log(0/H)^-gg is a bit lower than 
the other two, less than ~0.1dex lower. In this paper, 
we only plot the calibration results of 12±log(0/H)fl 23 
versus line ratios as representatives. 

The derived calibrations in this paper are only appro- 
priate for the metal-rich galaxies (12±log(0/H)>8.4), 
because the metallicities of these SDSS sample galax- 
ies are in the range of 8.4<12+log(0/H)<9.3. The cor- 
responding ranges of the various related strong-line ra- 
tios will be given in consequent subsections and Ta- 
bid. We omit the galaxies with lower metallicities, 
12±log(0/H)<8.4, in calibrations, and leave them to fu- 
ture work. 

3.1. Aperture effects in the SDSS 

Aperture bias is clearly important when comparing 
metallicities measured in the central regions of galax- 
ies to "global" quantities such as the luminosity (Kewley 
et al. 2005). The main question here is if the relation 
between the metallicity and the ionization parameter is 
the same in nuclear and global spectra. To estimate the 
aperture effect on the metallicity estimates, we compare 
our results about the SDSS galaxies with those of the 
local Nearby Field Galaxy Survey (NFGS) galaxies from 
Jansen et al. (2000a,b). Jansen et al. (2000b) have pub- 
lished the emission-line fluxes from both the integrated 
and nuclear spectra of their galaxies. We calculated their 
12±log(0/H) vs. [N 11] /Ha relations from their inte- 
grated and nuclear spectra, respectively, by using the 
same method as we adopted for the SDSS galaxies. 

We select 38 galaxies from the sample of Jansen et 
al. (2000b) following the criteria: 1) having the related 
fluxes for both the integrated and nuclear spectra avail- 
able; 2) confirming that they are star-forming galaxies by 
adopting the same criterion of Kauffman et al. (2003b) 
as we used for the SDSS galaxies; 3) using the same cal- 
culation method as for the SDSS galaxies, for example, 
for estimating Ay and deriving (O/H) abundances etc. 
Fig. El shows the comparison with the NFGS galaxies 
from the integrated spectra (the stars) and the nuclear 
spectra (the open squares). It is clear that all the data 
points from the nuclear spectra follow the SDSS galaxies 
very well, but some data points from the integrated spec- 
tra show lower 12±log(0/H) at a given [N il]/Ha value, 
~0.1-0.2dex lower. Fig.|5]may suggest a relatively lower 
ionization parameter for the outskirt regions of galaxies, 
but this conclusion is limited by the small number of 
galaxies with integrated spectra. Therefore, these cali- 
brations from SDSS galaxies may be more suitable to the 
nuclear spectra of galaxies. 

3.2. 1 2+log (O/H) versus [Nil] /Ha 
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The [N n]/Hct ratio is sometimes used to estimate 
metallicities of galaxies since it is not affected very much 
by dust extinction due to their close wavelength posi- 
tions (see Pettini & Pagel 2004). The near infrared spec- 
troscopic instruments, which were developed, can gather 
these two lines from the galaxies with intermediate and 
high redshifts. 

Metallicity calibrations of log([N n]A6583/Ha) (the 
N2 index) can be found in Storchi-Bergmann, Calzetti, 
& Kinney (1994) and Raimann et al. (2000). Re- 
cently, Denicolo et al. (2002) presented an analyt- 
ical linear formula by fitting the relation of a com- 
bined sample of 108 metal-poor galaxies (their O/H 
abundances were estimated from T e ) and 128 metal-rich 
galaxies (their O/H abundances were estimated from 
i?23 or S 23 (= ([S n]6717,6731+[S in]9069,9531)/H/3): 
12 + log (O/H) = 9.12 + 0.73 x N2. The data of Perez- 
Montero & Diaz (2005) are generally consistent with 
this calibration. Pettini & Pagel (2004) obtained a sim- 
ilar linear fit for a sample of 137 galaxies. Most of 
their samples have 7.0< 12+log(0/H)<8.6 based on T e 
measurements, and six (four are metal-rich) have oxy- 
gen abundances obtained from a detailed photoioniza- 
tion model. The linear fit given by Pettini & Pagel 
(2004) is: 12 + log (O/H) = 8.90 + 0.57 x N2, and an- 
other third-order polynomial fit with slight revision is: 
12+log(0/H) = 9.37+2.03xiV2+1.26x7V2 2 +0.32xiV2 3 
(valid in the range of -2.5 < N2 < -0.3). 

Here we present the relation of N2 index to (0/H)# 23 
for the 38,478 SDSS star- forming galaxies, and derive a 
corresponding calibration from these metal-rich galaxies. 
Fig-Ek shows our results for 12+log(0/H) j R 23 vs. N2 in- 
dex relations (the small blue points, without reddening- 
correction for [N n]/Ha ratios). The large red squares 
represent the 29 median N2 values in bins of 0.025 dex in 
12+log(0/H) within the range of 12+log(0/H)=8.5 to 
9.3. The 12+log(0/H) shows a clear increasing trend fol- 
lowing the increase of the N2 index, up to 12+log(0/H)~ 
9.0, and then, the galaxies having 12+log(0/H)> 9.0 
show a slightly decreasing N2 index. This trend can be 
understood as follows: when the secondary production of 
nitrogen dominates, at somewhat higher metallicity, the 
[N n]/Ha line ratio continues to increase, despite the de- 
creasing electron temperature; eventually, at still higher 
metallicities, nitrogen becomes the dominant coolant in 
the nebula, and the electron temperature falls sufficiently 
to ensure that the nitrogen line weakens with increasing 
metallicity (KD02). This is the first large sample to show 
this turnover relation between N2 and 12+log(0/H) at 
the very metal-rich region. An analytical third-order 
polynomial fit to the median-value points can be given 
in Fig. |2t as the solid line: 

12 + log(0/H) = aO + al x N2 + a2 x N2 2 

+ a3xJV2 3 . (3) 

The fitting coefficients a0, al, a2 and a3 are given in 
Table 1, as well the rms derivation (<r=0.058dex) for the 
data points from the fitting relation. This fit is obtained 
from the data points with -1.2 < N2 < -0.2 and 8.4< 
12+log(0/H)<9.3. The two dashed lines describe the 
discrepancy of 2a. A linear fit does not work well for 
these relations. 

However, we should notice that this calibration could 
result in large uncertainty when 12+log(O/H)>9.0 be- 



cause of the turnover of N2 index versus O /H in the very 
metal-rich region. Thus, we re-derive a third-order poly- 
nomial fit for this calibration for the part of sample galax- 
ies having 8.4<12+log(O/H)<9.0. We also re-derive the 
calibrations of the N2 index to 12+log(0/H)jf99 (the 
dot-long-dashed line in magenta color in Fig. 0t) and 
12+log(O/H)T04- All the fitting coefficients and the rms 
derivation a values are given in Table 1. 

In Fig. 13k we compare our calibration of [N n]/Ha vs. 
12+log(0/H) with other two calibrations published re- 
cently, as shown by the long-dashed and the dotted lines 
in Fig. Et- The long-dashed line represents the linear 
least squares fit given by Denicolo et al. (2002) for the 
whole sample of their 236 galaxies. The dotted line refers 
to the linear least squares fit obtained by Pettini & Pagel 
(2004) from their 137 sample galaxies. Both of these two 
calibrations show relatively lower O/H abundance com- 
pared to ours at a given N2 index, especially at the very 
metal-rich part. 

The difference from Pettini & Pagel (2004) is easily un- 
derstood since they only select metal-poor galaxies with 
T e measurements, except four metal-rich and two metal- 
poor ones. The fit for these metal-poor galaxies can be 
directly extrapolated to the metal-rich region and re- 
sults in a lower O/H at a given [N Il]/Ha value there. 
This may be related to the suggestion of Kennicutt et 
al. (2003), who pointed out that the abundance esti- 
mates from T e method are systematically lower by 0.2- 
0.5 dex than those derived from the strong-line "empiri- 
cal" abundance indicators, when the latter are calibrated 
with theoretical photoionization models. 

Denicolo et al. (2002) analyzed their metal-poor and 
metal-rich galaxies together and used the calibration of 
Kobulnicky et al. (1999) for converting i?23 to O/H for 
most of the metal- rich galaxies. Their fig.l shows that 
most of the metal-rich samples indeed lie on the left-top 
side of the fitting line. If only their metal-rich galaxies 
were considered to derive a calibration, it will predict a 
higher oxygen abundance at a given [N n]/Ha ratio than 
that they derived from the combined sample. To be clear, 
we re-derive the [N n]/Ha and 12+log(0/H) values for 
their metal-rich sample galaxies taken from Terlevich et 
al. (1991), in which the fluxes of the related lines were 
given as a table. The open triangles in Fig. present 
the H II galaxies from Terlevich et al. (1991) by using 
the same calibration as ours for the SDSS galaxies. The 
two samples show consistency. And the data of Terlevich 
et al. (1991) show more scatter, perhaps it comes from 
the different slit widths and resolutions used in different 
runs with different telescopes. 

3.3. 12+log(0/H) versus [O m]/[N n] 

Alloin et al. (1979) firstly introduced the quantity 
03N2 = log{([0 in] A5007/H/3)/([N n] A6583/Ha)} 5 . 
Recently, Pettini & Pagel (2004) presented one such 
calibration by doing linear least squares fitting for a 
sample of 65 (from the 137) galaxies in the range of 
-1 < 03^2 < 1.9. They obtained a calibration of 
12 + log (O/H) = 8.73 - 0.32 x 03N2. 

We can obtain such a calibration from the large sam- 

5 This definition is slightly different from the original one pro- 
posed by Alloin et al. (1979) who included both [O III] doublet 
lines in the numerator of the first ratio. 
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pie of SDSS galaxies. Fig. |2t> shows the i?23-derived 
12+log(0/H) versus 03N2 relations for our sample 
galaxies. The solid line gives a linear least squares fit: 

12 + log (O/H) = 60 + 61 x 03N2, (4) 

with the coefficients bO and bl given in Table 1, as well as 
the rms derivation. The two dashed lines show the dis- 
crepancy of 2cr. The dotted line refers to the fit of Pettini 
& Pagel (2004) for their 65 metal-poor galaxies, which 
gives lower O/H than that of our samples at a given 
03N2 value; about 0.1 dex lower. The reason may be 
that the galaxies with T e -abundances were mainly used 
in their study. 

To get a better fit for the details of the relation, a third- 
order polynomial fit is obtained and shown as the long- 
dashed line in Fig. [3Jd, which revises the linear equation 
slightly. The corresponding fitting coefficients of aO, al, 
a2 and a3 are given in Table 1 with the same meaning 
as Eq.(l), as well as the rms derivation. These linear 
and three-order polynomical calibrations were obtained 
from the sample galaxies with 8.4< 12+log(0/H) <9.3 
and -0.7 <03N2< 1.6. The obvious advantages of using 
the [N n] and [O in] lines is that they are unaffected by 
absorption lines originating from the underlying stellar 
populations, they lie close to Balmer lines that can be 
used to eliminate errors due to dust reddening, and they 
both are strong and easily observable in the optical band. 

3.4. 12+log(0/H) versus [Nu]/[0 n] 

We can derive the calibration of [N il]/[0 n] ratio to 
12+log(0/H) from the large sample of the SDSS galax- 
ies. Fig. yfc presents these data points in the relation of 
12+log(0/H) versus [N ll]/[0 n]. A linear least squares 
fit reproduces the data distribution well, as shown by the 
solid line. The two dashed lines show 2a discrepancies. 
A third-order polynomial fit will revise the fitting equa- 
tion slightly. The coefficients of the two fits and the rms 
derivation a of the data to the line are given in Table 1, 
as well as the appropriate ranges of 12+log(0/H) and 
log([N n]/[0 ii]). The fitting formulas follow the same 
format as Eq.(2) and Eq.(l). 

Fig. Et shows that the [N il]/[0 n] ratios obviously 
increase with the increasing metallicities, which can be 
explained by the predominant secondary enrichment of 
[N n] above metallicities of 12+log(0/H)>8.6 (Alloin ct 
al. 1979; KD02), and the stronger decrease in the number 
of collisional excitations of the blue [O n] lines relative 
to the lower energy [N n] lines due to the low electron 
temperature in high metallicity environments. 

3.5. 12+log(0/H) versus [Nu]/[S Ii] 

We can derive the calibration of [N n]6583/[S n]6717, 
6731 to 12+log(0/H) by using the large sample of SDSS 
galaxies. Fig. [Sji shows these observed results, which can 
be explained by a linear least squares fit and is given as 
the solid line. The two dashed lines show the 2a discrep- 
ancy. To get a more detailed fit, a third-order polynomial 
fit has been obtained (the long-dashed line in Fig. |3Ji), 
which revises the relation slightly. The coefficients of 
the two fits and the rms derivations are given in Ta- 
ble 1, as well the appropriate ranges of 12+log(0/H) and 
log([N ii] / [S ii]). The fitting formulas follow the same 
format as Eq.(2) and Eq.(l). 



It shows that the 12 + log(0/H) abundances of the 
sample galaxies increase with their log([N n]/[S n]) ra- 
tios. The reason may be that at high metallicity, nitrogen 
is a secondary nucleosynthesis element and sulphur is a 
primary process element (KD02). Both lines have simi- 
lar excitation potential since they are close to each other 
in wavelength. 

3.6. 12+log(0/H) versus [S h]/Hq 

Although [S n]/Ha ratio can be easily derived from 
[N Il]/Ha and [N n] / [S ii] ratios, it is still interesting to 
directly present the [S n]/Ha to 12+log(0/H) relations 
for this large sample of SDSS galaxies. Fig. ^ shows 
that these data points show a double- valued distribution: 
because of the turnover around 12+log(0/H)^8. 9-9.0, 
there will be double- valued oxygen abundances at a given 
[S n]/Ha ratio. Thus, it is difficult to fit a simple and 
useful analytic formula from these observational data. 

We should notice that the [S ll]/Ha ratio has never 
been suggested or used as a metallicity indicator, because 
it is far more sensitive to ionization (see Sect. 4) than to 
metallicity, and it results in double- valued (O/H). If it is 
able to measure the [S n] and Ha lines, it is likely that 
we have also obtained [N ii], and then we could use the 
[N n]/Ha calibration. 

3.7. 12+log(0/H) versus [O lll]/H/3 

The [O m]4959,5007/H/3 ratios can also be used to 
estimate the metallicities of galaxies. This was called 
"i?3 method" (Edmunds & Pagel 1984), and was shown 
as: R 3 = 1.35 x (I[o///] 5007 and then, log(0/H) = 
-0.69 x logi? 3 - 3.24, with a range of -0.6 < logi? 3 < 1.0 
(Vacca & Conti 1992). 

Here we present a new calibration for this ratio vs. 
O/H derived from the SDSS sample galaxies (Fig. Ef). 
The linear least squares and the third-order polynomial 
fits are obtained, and given as the solid and the long- 
dashed lines in Fig. [3f, respectively. The coefficients of 
the two fits and the rms derivations are given in Ta- 
ble 1 following Eq.(2) and Eq.(l), as well the appropriate 
ranges of 12+log(0/H) and logi?3. The two dashed lines 
show the 2a discrepancy. 

The observated results of this large sample of 
SDSS star-forming galaxies and the derived ana- 
lytic calibrations show that the higher [N ll]/Ha, 
[N n]/[0 ii] and [N n]/[S ii], and the lower [O m]/H/3, 
([O m]/H/3)/([N ll]/Ha) ratios will result in the higher 
12+log(0/H) abundances generally. The rms derivations 
of the observational data relative to the calibration re- 
lations are in 0.032-0.078 dex. [N n]/[0 ii] has the low- 
est rms derivation. These calibrations from strong-line 
ratios to metallicities can be used to calibrate oxygen 
abundances of galaxies in future studies. 

4. COMPARISON WITH THE PHOTOIONIZATION MODELS 

Our strong-line calibrations are in good qualita- 
tive agreement with those of KD02. For example, 
both calibrations show a turnover of [N n]/Ha near 
12+log(O/H)~9.0. One major difference is that KD02 
provided strong-line calibrations for seven different val- 
ues of the ionization parameter, whereas we assume no 
explicit ionization parameter dependence. Our approach 
will be less accurate in some cases, but it is necessary 
when only a few emission lines can be measured and there 
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is too little information to determine both ionization pa- 
rameter and metallicity. It is also unclear that the full 
range of ionization parameters used in the KD02 mod- 
els are required by the data. We explore this in Fig. El 
where we plot the relations between R23 and other strong 
line ratios in the SDSS data and overlay the KD02 model 
grids. Fig. also gives the similar information. 

In Figs. Efc-f, the y-axis log(i?23) is plotted from the 
lower to the higher values, which corresponds to an 
increasing trend of metallicity for the galaxies in the 
metal-rich branch. The seven lines in each plot repre- 
sent the model results of KD02 for ionization parameters 
g = 5 x 10 6 , 1 x 10 7 , 2 x 10 7 , 4 x 10 7 , 8 x 10 7 , 1.5 x 10 8 , 
and 3 x 10 8 cms -1 , respectively. From the dotted line to 
the solid line, the ionization parameters increase orderly. 
The lines in Figs. EM start from 12+log(0/H)=8.4 for 
the metal-rich branch, and the turn-over trend of the 
model grids at the low-metallicity end corresponds to 
the turnover of log(i?23) itself from the metal-rich to the 
metal-poor branches. 

- [N 11] /Ha Fig. Et shows the comparison for 
log(i?23) versus log([N n]/Ha). It shows that these ob- 
served points are reasonably consistent with the model 
results of KD02. The increasing trend of [N 11] /Ha fol- 
low the decreasing log(i?23)a-nd then the turnover from 
about log(i?23)~ 0.5 correspond to the increasing and 
then the turnover of [N 11] /Ha following the increasing 
12+log(0/H) abundances as shown in Fig.Et. However, 
[N II] /Ha is sensitive to the ionization parameter. It 
seems that q = 1.5 X 10 8 and 3 x 10 8 cms -1 are too 
high and q — 5 x 10 6 cms -1 is too low to explain the 
distribution of these star-forming galaxies. Most of the 
sample galaxies fall in the range of ionization parameter 
q from 1 x 10 7 to 8 x 10 7 cms -1 , which is narrower than 
the range proposed by the models. 

[N 11] /[O ill] Fig. shows the compari- 

son between the observed relations of log(i?23) versus 
log([N 11] /[O ill]) with the photoionization model results 
of KD02. It shows that the observed relations are con- 
sistent with the model results: the log([N 11] /[O ill]) ra- 
tios of these galaxies increase following the decreasing 
log(i?23) (the increasing metallicities). It also shows that, 
most of the data points distribute in a range of ionization 
parameter q from q = 1 x 10 7 to 8 x 10 7 cms -1 , which 
is similar to log([N Il]/Ha). [N 11] /[O ill] is sensitive to 
the ionization parameter. 

-[N II] /[O 11] Fig. Et shows the comparison 

between the observed results and the photoionization 
model results of KD02 for the relations of log(i? 2 3) ver- 
sus log([N 11] /[O 11]). It shows that the [N 11] /[O 11] 
ratios of the sample galaxies increase following the de- 
creasing log(i?23) (the increasing metallicities), which is 
very consistent with the general trend of the model re- 
sults. As KD02 mentioned, [N 11] /[O 11] is almost in- 
dependent of the ionization parameter. Our large set 
of data points shows a very narrow distribution in the 
log(i?23) vs. [N n]/[0 11] relation, which is very consis- 
tent with KD02's model results. However, [N n]/[0 11] 
ratios are affected strongly by dust extinction inside the 
galaxies. We will discuss this in Fig. |S| and in the last 
paragraph of this section. 

[N 11] /[S 11] Fig. EJi shows the comparison 

between the observed results and the photoionization 
model results of KD02 for the relations of log(i?23) ver- 



sus log([N n]/[S 11]). It shows that the [N n]/[S 11] ra- 
tios of the sample galaxies increase with the decreasing 
log(i?23) (the increasing metallicities), which is consis- 
tent with the general trend of the models. However, the 
fact is that these SDSS data points show relatively lower 
log(i?23) (higher metallicities) than the model predicts at 
a given [N n]/[S 11] ratio, about 0.2 dex discrepancy. One 
of the reasons may be that the amount of S going into 
[S in] compared to [S 11] is not well understood, so this 
could introduce some uncertainties into the photoioniza- 
tion models (see KD02). KD02 commented that their 
[N Il]/[S II] technique will systematically underestimate 
the abundance by ~0.2dex compared to the average of 
the comparison techniques. The comparison of models 
and data shows that the models fail to reproduce the 
strengths of the [S 11] lines, so we cannot use the mod- 
els to infer abundances or ionization parameters from 
[N n]/[S 11], as well [S 11] /Ha (see below). These high- 
light the need for empirical calibration. 

- [S 11] /Ha Fig. Et shows the comparison be- 
tween the observed results and the photoionization model 
results of KD02 for the relations of log(i?23) versus 
log([S 11] /Ha). There is a weak correlation between 
these two parameters, meaning that log([S 11] /Ha) de- 
creases following the decreasing log(i?23) with a very 
weak turnover around log(i?23)~0.5. For the [S n]/Ha 
calibration, the SDSS data points agree with a relatively 
higher model q value than in the [N 11] /[S 11] case. 

[O iii]/H/3 Fig. Et shows the comparison 

between the observed results and the photoionization 
model results of KD02 for the relations of log(i?23) 
versus log([0 iii]/H/3). The model results of [O in]/H/3 
from KD02 were deduced using their three other ratio 
parameters, i.e. log([N 11] /[O in]), log([N u}/[0 11]) 
and log(i?23), since KD02 did not give a direct formula 
for this calibration. These SDSS sample galaxies 
follow the trend of the model results well, meaning the 
log([0 in]4959,5007/H/3) ratios decrease following the 
decreasing log(i?23) (the increasing metallicities). This 
ratio depends strongly on the ionization parameter. 
Most of the data points are in a range of q = 1 X 10 7 
cms -1 to 8 x 10 7 cms -1 , except in the very metal-rich 
region with log(i?23)<0.2. 

Figs.EKjbif show that these star- forming galaxies have 
a relatively narrower range of ionization parameter q 
than the range proposed by the photoionization mod- 
els of KD02. The actual range of q of these galaxies 
is 1 x 10 7 — 8 x 10 7 cms -1 , narrower than the range of 
models, q = 5 x 10 6 — 3 x 10 8 cms -1 . We do not con- 
sider Figs. Et,d,e for this range since [N n]/[0 11] is in- 
sensitive to q and the models provide a poor match to 
[N n]/[S 11] and [N 11] /Ha (KD02). To show more clearly 
the range of q of these sample galaxies, we compare their 
observed results with the model grids in Fig. OD for the 
[O m]/[N 11] versus [N 11] /[O 11] relations following Do- 
pita et al. (2000) (their fig. 7). The model grids are taken 
from Kewley et al. (2001) 6 , their instantaneous zero-age 
starburst models based on the PEGASE spectral energy 
distribution (SED), which are the most appropriate to 
use to model H 11 regions (KD02), and show informa- 

6 Thanks to Lisa Kewley for putting their model grids on the 
web: http://www.ifa.hawaii.edu/~kewley/Mappings/ 



Oxygen abundance calibrations and N/O abundances of SDSS galaxies 



7 



tion about both metallicity and q. This diagram shows 
that these nearby star-forming galaxies have metallicities 
Z = 0.2 — 2.0Zq, and the actual ionization level of their 
ionized gas have ionization parameters in the range of 
q = 1 x 10 7 to 8 x 10 7 cms _1 with some exceptions at the 
very metal-rich end. 

Indeed, the parameter q is related to the ionization 
parameter U by dividing the speed of light c: U = q/c, 
which is more commonly used (e.g. McGaugh 1991). 
Therefore, the ionization parameters of KD02's models 
correspond to a range from logJ7=— 3.78 to logJ7=— 2. 
Our Figs.0K,b,f show that the actual range of ionization 
parameters of these SDSS sample galaxies are from about 
\ogU=— 3.5 to —2.5. This range is also consistent with 
what found by Dopita et al. (2000) for the extragalactic 
H II regions. 

We also compare these SDSS star-forming galaxies 
with the dwarfs from van Zee & Haynes (2006, the filled 
triangles) and local spirals from van Zee et al. (1998, 
the open circles) in Fig. on the basis of the relation of 
[N n]/[0 n] versus logi? 2 3- Although [N n]/[0 n] is in- 
sensitive to ionization parameters, log7?23 is, so this plot 
can help to understand the evolutionary status of these 
sample galaxies. It is clear that these SDSS galaxies 
follow the metal-rich spirals well, but the dwarfs probe 
a metallicity range not covered by the SDSS galaxies- 
only very few dwarfs have 12+log(0/H)>8.4. 

We should notice that dust extinction will seriously 
affect the metallicity calibration from the line ratio of 
[N n]/[0 II] since the blue line [O II] is strongly af- 
fected by dust extinction. Fig. [7| shows the results with- 
out extinction correction for [N n]/[0 n] for our sample 
galaxies, which shows a much wider distribution, and 
is quite far moved from KD02's model results. The 
comparison between Fig. [7| and Fig. QJj shows that: 1) 
the dust extinction in these SDSS galaxies could not 
be omitted; 2) the extinction coefficients of the individ- 
ual galaxies are quite different, thus, the data can show 
large scatter when the dust extinction inside the individ- 
ual galaxies were not estimated properly. For other ra- 
tios, like ([O m]/H/3)/([N n]/Ha), [N n]/Ha, [N n]/[S n] 
and[0 iii]/H/3, dust extinction does not have much effect 
on them due to the close positions of the two related lines 
in wavelength. We did not consider the dust extinction 
for [N il]/Ha and [N n]/[S II ] ratios in our calculations. 

5. NITROGEN-TO-OXYGEN ABUNDANCE RATIOS 

It is possible to estimate the N abundances of galax- 
ies from strong optical emission lines, which can help to 
understand the origin of nitrogen. 

The basic nuclear mechanism to produce nitrogen is 
well understood - it must result from CNO processing of 
oxygen and carbon in hydrogen burning. However, the 
nucleosynthetic origin of nitrogen has a "primary" and a 
"secondary" component, which is still in debate. If the 
"seed" oxygen and carbon are those incorporated into 
a star at its formation and a constant mass fraction is 
processed, then the amount of nitrogen produced is pro- 
portional to the initial heavy-element abundance, and 
the nitrogen synthesis is said to be "secondary". If the 
oxygen and carbon are produced in the star prior to the 
CNO cycling (e.g. by helium burning in a core, followed 
by CNO cycling of this material mixed into a hydrogen- 
burning shell), then the amount of nitrogen produced 



may be fairly independent of the initial heavy-element 
abundance of the star, and the synthesis is said to be 
"primary" (Vila-Costas & Edmunds 1993). In general, 
primary nitrogen production is independent of metallic- 
ity, while secondary production is a linear function of 
it. 

From a theoretical point of view, the secondary pro- 
duction of nitrogen should be common to stars of all 
masses, whereas the primary production should arise 
only from intermediate-mass stars (4 < M/Mq < 8) un- 
dergoing dredge-up episodes during the asymptotic giant 
branch evolution. In particular, when the third dredge- 
up is operating in conjunction with the burning at the 
base of a convective envelope (hot-bottom burning), pri- 
mary nitrogen can originate (Renzini & Voli 1981; Mat- 
teucci 1986). The large set of database of galaxies from 
SDSS must provide important information on the nucle- 
osynthetic origin of nitrogen. The N/O ratio as a func- 
tion of O/H is the basic method to study the N abun- 
dance of galaxies. 

5.1. The method and results 

Firstly, we use the formula given by Thurston et al. 
(1996) to estimate the electron temperature in the [N n] 
emission region (T[nh]) from logi?23: 

T [NII] = 6065 + 1600a; + 1878a; 2 + 2803a; 3 , (5) 

where a;=log(R23), TJ N n] is in units of K. This for- 
mula was derived from the data points with about - 
0.6<logi?23<0.8, and 5000K<T [NII] <10000K (see fig.l of 
Thurston et al. 1996). 

Then, log(N/0) values are estimated from the ([N n] 
AA6548, 6583)/ ([O n]A3727) emission-line ratio and T [NII] 

temperature by assuming ^ = and using the con- 
venient formula based upon a five-level atom calculation 
given by Pagel et al. (1992) and Thurston et al. (1996): 

M£) = log [NI S 8 ra 6 f 3 + 0.307 - 0.0210^,,, 
0.726 



where twii] = 10 4 7] NII ]- We consider the flux of 
[N n]A6548 is equal to one-third (0.333) of that of the 
[N li]A6583 in calculations. 

Fig. H gives the log(N/0) vs. 12+log(0/H) relations 
for our sample galaxies. We also plot models for the 
primary and secondary origin of nitrogen taken from 
Vila-Costas & Edmunds (1993): the dot-dashed line 
refers to the "primary" component, the long-dashed line 
refers to the "secondary" component, and the solid line 
refers to the combination of the two components. The 
log (N/O) abundances of the sample galaxies are consis- 
tent with the combination of the "primary" and "sec- 
ondary" components, but the secondary one dominates 
in these metal-rich galaxies. This result is similar to 
the previous studies, e.g. Shields et al. (1991), Vila- 
Costas & Edmunds (1993), Contini et al. (2002), and 
Kennicutt et al. (2003). In addition, 67 H n regions 
in 21 dwarf irregular galaxies taken from van Zee &; 
Haynes (2006) are also plotted in Fig.|Hlas triangles. The 
log (N/O) abundances of these dwarf irregulars are dom- 
inated by the "primary" component. The SDSS galaxies 
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with 8.0<12+log(O/H)<8.5 in the metallicity transition 
region are also presented as the stars in Fig. El They are 
just distributed in the region between the H n regions in 
dwarf irregulars and the metal-rich galaxies. 

5.2. The analyses and discussions 

The SDSS sample galaxies show increasing log(N/0) 
abundance ratios following increasing 12+log(0/H) 
abundances. The behavior of N/O with increasing metal- 
licity provides clues about the chemical evolution history 
of the galaxies and the stellar populations responsible 
for producing oxygen and nitrogen. Oxygen is mainly 
produced in short-lived, massive stars, and is ejected 
to the ISM by Type II supernova explosions (SNe). In 
addition, nitrogen is mainly produced in the long-lived, 
intermediate- and low-mass stars, and is ejected to ISM 
through stellar wind. In particular, the burning at the 
base of the convective envelope (hot-bottom burning) 
and in conjunction with the third dredge-up process in 
intermediate-mass stars contribute an important fraction 
of nitrogen (Renzini & Voli 1981; van den Hoek & Groe- 
newegen 1997; Liang et al. 2001; Henry et al. 2000). 
Therefore, nitrogen will be delayed in its release into the 
ISM as compared to oxygen. 

Contini et al. (2002) analyzed the delayed release of ni- 
trogen. As they discussed, during a long period of quies- 
cence, intermediate-mass star evolution will significantly 
enrich the galaxy in nitrogen but not in oxygen; during 
starburst, however, N/O drops while O/H increases as 
the most massive stars begin to die and supernovae re- 
lease oxygen into the ISM. A few tens of Myr after the 
burst, the massive stars producing oxygen will be gone. 
Then N/O will rise again as intermediate-mass stars be- 
gin to contribute to the primary and secondary produc- 
tion of nitrogen. A galaxy undergoing successive star- 
bursts will make the N/O ratio increase following the in- 
creasing O/H abundance. However, the increasing N/O 
following O/H in the metal-rich region does not have to 
require star formation to occur in periodic starbursts. 
Rather, as long as the recent star formation activity is 
less than the past star formation rate, the delayed release 
of nitrogen from the aggregate intermediate mass stellar 
population will slowly increase the N/O ratio since it is 
not balanced by additional production of oxygen in high 
mass stars (van Zee & Haynes 2006). 

One fact in Fig. [S] should be noticed: when the 
O/H increases to be 12+log(O/H)~9.0, the slope of 
the increase of log(N/0) becomes steeper, i.e. the 
log(N/0) shows an upturn trend at higher metallicities 
than 12+log(O/H)~9.0, or O/H then saturates for the 
increasing N/O. The main reason for this may be that the 
sensitivity of R23 to metallicity saturates somewhat at a 
high metallicity. To understand more about this, we red- 
erive the relations of log(N/0) versus the other two sets 
of i?23-derived oxygen abundances, 12+log(0/H)^-gg and 
12+log(0/H)z94. The later one is derived from the con- 
version of R23 to 12+log(0/H) given by Zaritsky et al. 
(1994). They show the same trend as ^-Mog^/H)^. 
We then take the Bayesian estimates of 12+log(O/H)T04 
from Tremonti et al. (2004) to replace the i?23-derived 
oxygen abundances to present the N/O vs. O/H rela- 
tions of the sample galaxies. The results are given in 
Fig. El The saturation of 12+log(0/H) to log(N/0) be- 
comes much weaker now. This confirms that the satura- 



tion of O/H to N/O at the metal-rich end mainly comes 
from the insensitivity of R23 to 12+log(0/H) there. In- 
deed, our Fig^has shown this saturation of i?23-derived 
oxygen abundances at the high metallicity end. That 
may mean that the upturn in N/O at high abundance is 
not physical. 

From the comparisons above, we should note that care 
should be taken in interpreting features in the N/O 
vs. O/H diagram because of systematic uncertainties in 
metallicity calibration. If we use T e rather than strong- 
line calibrations, O/H would be lower but N/O would 
not be much affected since it has only a weak depen- 
dence on T e . The n et ef fect would be a shift of the data 
to the left in Figs. 18191 The data points will then be 
more consistent with the H 11 regions in M101 studied by 
Kennicutt et al. (2003) in the relations of N/O vs. O/H 
diagram. The magenta dotted line in Fig.|S]is the "best" 
model B of the numerical one-zone models for the evo- 
lution of N/O versus O/H given by Henry et al. (2000). 
Kennicutt et al. (2003) found that the predicted N/O 
mode by this model is too low at intermediate abun- 
dance (12+log(0/H)T e =8.2-8.6) to predict the observed 
results of the H 11 regions in M101. Our SDSS sam- 
ple galaxies show a similar discrepancy from the model: 
at the intermediate abundance 12+log(O/H)<9.0, the 
model predicts too low N/O at a given O/H; for the 
12+log(O/H)>9.0, although the model predicts an in- 
creasing trend of N/O following O/H, the predicted N/O 
is still generally lower than the observed results at a given 
O/H. Indeed, Henry et al. (2000) themselves already 
pointed out that their model B provides a fit to the NO 
envelope compared with the observational data (see their 
fig.3b). 

In addition, dust extinction is a very important factor 
affecting the estimated (N/O) abundances, as well as the 
O/H abundances. It can be understood from Fig. (by 
comparing with Fig. Ufc). The dust extinction of these 
SDSS sample galaxies has been considered when deriving 
their abundances. Moreover, it is a homogeneous sam- 
ple taken from one single facility using the same set-up, 
which minimizes the systematic errors from observation 
runs, and the [N 11] /[O 11] ratio is almost independent of 
the ionization parameter (KD02). Therefore, these SDSS 
data points would not show much spread in the N/O 
vs O/H relations due to observational errors. The ob- 
servational scatter of these galaxies are about 0.062 dex 
estimated from Fig. ED 

The scatter in the observed log(N/0) vs. log(0/H) 
relation is relatively small. Nevertheless, we can use our 
large dataset to see if the residuals correlate with any 
physical properties of galaxies. If the release of nitro- 
gen is indeed delayed with respect to oxygen, we might 
expect to find a correlation between the N/O residu- 
als and the stellar birthrate, which is the ratio of the 
present to past average star formation rate (SFR). The 
equivalent widths of Ha emission lines, log(EW(Ha)), 
provides a good proxy for the stellar birthrate (Kenni- 
cutt et al. 1994). Fig. 1101 shows the relation between 
log (N/O) residual and the log(EW(Ha)) for the sample 
galaxies. log(N/0) res iduai is the N/O abundance hav- 
ing its dependence on O/H abundances been removed. 
To do so, first we obtain a linear least squares fit for 
the relation between log(N/0) vs. 12+log(0/H)ro4 for 
the sample galaxies, then the observed log(N/0) abun- 
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dances were subtracted by the fitted values predicted by 
this linear relation. Fig. ED shows that the log(N/0) 
abundances decrease following the increasing EWs(Ha) 
for these galaxies. The solid line refers to a linear least 
squares fit for this observed relation with a rms ^0.057: 
N 

log(Q)rcsiduai = -0.102 x log(EW(Ha)) + 0.140. (7) 

This relation means that the log(N/0) increases follow- 
ing the decreasing star formation rates, which is consis- 
tent with the suggestion of van Zee & Haynes (2006) for 
the relation between log(N/0) and B-V colors for their 
dwarf galaxies. They suggested that the high N/O ratio 
are correlated with redder systems and decreasing star 
formation rates. 

6. CONCLUSION 

We have use the observational data of a large sample 
of 38,478 star-forming galaxies selected from SDSS-DR2 
to derive oxygen abundance calibrations of 12+log(0/H) 
versus several metallicity-sensitive emission-line ratios, 
including [N n]/Ha, [O m]/[N n], [N n]/[0 n], 
[N n] /[S n], [S n] /Ha, and [O m]/H/3. These calibra- 
tions are very useful when the "T e -method", the most 
accurate method, and the empirical "i?23-niethod" are 
not available to estimate oxygen abundances of galax- 
ies. In addition, they can overcome the "double-valued" 
drawback of the i?23-method. 

We fit the observed relations with a linear least squares 
fit and/or a third-order polynomial fit. All the fitting co- 
efficients and the rms derivations (0.032-0.078 dex) have 
been given in Table 1. When we only select the sample 
galaxies with higher signal-to-noise ratio emission lines, 
e.g. >10er, the rms derivations only have small changes, 
which means that probably most of the scatter is intrinsic 
and not an artifact of observational errors. The scatter- 
ing of the data points mainly come from the different 
ionization parameters in the individual galaxies, which 
is confirmed by the smallest rms value of [N n]/[0 n] ra- 
tio, which is independent on the ionization parameters. 
Such a large sample of galaxies can show some infor- 
mation never before obtained from much smaller sam- 
ples, for example, the turnover of the relations between 
12+log(0/H) versus [N n]/Ha at 12+log(O/H)>9.0. 

Among these calibrations, the [N n]/[0 n] ratio should 
be the best one for metallicity calibration in the metal- 
rich branch (same as KD02) because it shows a monoton- 
ical increase following the increasing metallicity and less 
scatter than other line ratios, and it is independent of the 
ionization parameter. However, dust extinction must be 
estimated properly before using this calibration indica- 
tor because of the blue line [O II] and the far positions 
of the two lines in wavelength. ([O m]/H/3)/([N n]/Ha) 
is also a good indicator, a bit better than [N n]/[S II] 
and [O iii]/H/3 due to the relatively smaller rms deriva- 
tion. The discrepancy between the data and models for 
[N n]/[S n] highlights the need for empirical calibrations. 
[N n]/Ha is only a good indicator for the galaxies with 
12+log(O/H)<9.0, and the problem at higher metallici- 
ties is the turnover of [N n]/Ha following increasing O/H. 
However, [N n]/Ha may be one of the most useful cal- 
ibrations since it is not affected by dust extinction due 
to the very close positions of the two related lines in 
wavelength, and can be detected by the near infrared in- 
struments for the intermediate- and high-z star-forming 



galaxies. We should notice that the derived calibration 
will result in a higher 12+log(0/H) abundance at a given 
[N n]/Ha ratio than those of Pettini & Pagel (2004) and 
Denicolo et al. (2002). The main reason for this differ- 
ence may be due to the discrepancy between the metal- 
licity estimates from the T e -method and i?23-method. 

The resulting calibrations are more suitable for lumi- 
nous and metal-rich galaxies since they have metallicities 
in the range 8.4< 12+log(0/H)fl 23 <9.3. Because our 
calibrations are derived from SDSS data, they will be 
best suited for data which samples the inner few kilopar- 
secs of galaxies rather than the whole disk. This aperture 
effect is most significant for line ratios with a strong de- 
pendence on ionization parameter (e.g. [N n]/[0 in]). 

The observed calibration relations are compared with 
the photoionization models from KD02, which gener- 
ally show good consistent trends, including the turnover 
trend of [N il]/Ha at the metal-rich end, the indepen- 
dence of [N n]/[0 Ii] on ionization parameter, etc. How- 
ever, most of these sample galaxies are distributed in a 
narrower range of the ionization parameter q than the 
models. The actual range of q is from 1 x 10 7 to 8 x 10 7 
cms^ 1 , i.e., log[/=-3.5 to -2.5, but the models cover from 
q = 5 x 10 6 to 3 x 10 8 cms" 1 , i.e., log[/=-3.78 to -2. 

Another contribution of this work is that we ob- 
tained the log(N/0) abundance ratios for this large 
sample of SDSS star-forming galaxies with 8.4< 
12+log(0/H)<9.3. Their log(N/0) abundances are con- 
sistent with the combination of the "primary" and "sec- 
ondary" components of nitrogen, but the "secondary" 
one dominates at high metallicities. The increasing 
log(N/0) abundances that follow the increasing metal- 
licities for these metal-rich galaxies can be the prediction 
of continuous but declining star formation rates, which is 
confirmed by the linear relation with a negative slope of 
the log(N/0)rcsiduai versus log(EW(Ha)). The scatter of 
these observational data are ^0.062 dex in the log(N/0) 
versus 12+log(0/H) relations. However, we should no- 
tice that dust extinction strongly affects the estimates 
of the log(N/0) abundance ratios and the corresponding 
scatter for the sample galaxies. 

In summary, this work provides several useful metal- 
licity calibrations based on strong emission-line ratios of 
38,478 SDSS star-forming galaxies, and estimates the 
log(N/0) abundance ratios for such a large number of 
galaxies. These can be used as references for future stud- 
ies of metallicities of galaxies. We should notice that, 
when these calibrations are used for high redshift ob- 
jects, the indices based on [N n] may have some extra 
uncertainties since for very high redshifts N/O vs. O/H 
may be different from the local case because the galax- 
ies are not old enough to produce significant amounts of 
secondary nitrogen. 
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Fig. 1. — Comparison between our metallicity estimates from the extinction-corrected R23 parameters (R23) and those Bayesian estimates 
from Tremonti et al. (2004) (T04) derived using the model of Chariot et al. (2006). These two estimates are very consistent. The solid 
line refers to the line of equality, and the ±0.15 dex discrepancies from the solid line are also given as the two dashed lines. 
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Fig. 2. — Comparison between our sample galaxies (the small blue points) and the NFGS galaxies at local from Jansen et al. (2000a, b): 
the stars refer to the results from the integrated spectra, and the open squares refer to the results from the nuclear spectra (see text). 
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Fig. 3. — Analytic calibration relations between 12+log(0/H) and several metallicity-sensitive strong-line ratios based on the observed 
relations of the SDSS sample galaxies (the small blue points): (a.) versus N2 index (without reddening correction for [Nn]/Ha). The 
large squares represent the 29 median values in bins of 0.025 dex in 12+ log(0/H) within 12+log(0/H) from 8.5 to 9.3. The open triangles 
refer the local H II galaxies from Denicolo et al. (2002) (from Terlevich et al. 1991). The dot-long-dashed line represents the calibration 
with 12+log(0/H)x99- (b.) versus 03N2 index relations, (c.) versus log([N ll]/[0 II]). (d.) versus log([N II] / [S II]) (without reddening 
correction for [Nll]/[S o]). (e.) versus log([S Il]/Ha). (f.) versus log([0 iii]/H/3) ratios. The solid lines refer to the fits to the data (see text), 
and the two dashed-lines show the discrepancy of 2ct; the long-dashed lines in Figs.b,c,d,f refer to the three-order polynomial fits for the 
observational relations. The dotted lines in Figs. a, b refer to the linear least-squares fit given by Pettini & Pagel (2004), the long-dashed 
line in Fig. a refers to the result of Denicolo et al. (2002). The dot-long-dashed line in Fig.f refers to the linear least-squares fit given by 
Vacca & Conti (1992). All the fitting coefficients for the SDSS galaxies and the rms derivations have been given in Table 1 following the 
format of Eq.(2) and Eq.(l). 
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Fig. 4. — Comparison between our SDSS sample galaxies {the small blue points) and the photoionization model results of KD02 
{the lines) for the relations of log(i?23) versus several strong-line ratios: (a.) log([N ii]/Hq); (b.) log([N ll]/[0 III]); (c.) log([N ll]/[0 n]); 
(d.) log([N II] /[S II]); (e.) log([S n]/Ha); (f.) log([0 m]/H/3). The lines represent the model results of KD02 with q = 5 X 10 6 , 1 X 10 7 , 
2 x 10 7 , 4 x 10 7 , 8 x 10 7 , 1.5 x 10 s , and 3 x 10 8 cm s , shown as from the red dotted line to the black solid line, respectively. The seven 
lines in Fig.f were deduced from other three line ratios, i.e. log([N ll]/[0 in]), log([N n]/[0 ll]) and log(i?23), since there are no direct 
log([0 iii]/H/3) calibration relation given in KD02. 




Fig. 5. — Comparison between the observed results of the SDSS sample galaxies and the photoionization models of Kewley et al. (2001) 
following fig. 7 of Dopita et al. (2000). The solid lines refer to constant ionization parameters, with the grids of q = 5 X 10 6 , 1 X 10 7 , 2 X 10 7 , 
4 X 10 7 , 8 X 10 7 , 1.5 X 10 8 , 3 X 10 8 cms -1 ; and the dotted lines refer to constant metallicity, with the grids of Z =0.1, 0.2, 0.5, 1.0, 1.5 
and 2.0Z Q . 
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Fig. 6. — Comparison between our sample galaxies and the spiral galaxies form van Zee et al. (1998) and the dwarf irregulars from van 
Zee & Haynes (2006). 




Fig. 7. — Comparison the log(i?23) versus log([N ll]/[0 II]) relations between our sample galaxies without correcting dust extinction for 
[N ll]/[0 II] and the model results of KD02. The comparison of this plot with Fig. |3Ji shows that the dust extinction strongly affect the 
[N n]/[0 II] ratios and their scatter, hence the calibrated metallicities. 
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Fig. 8. — log(N/0) abundance ratios of our SDSS sample galaxies as functions of their 12+log(0/H) abundances derived from R23 by 
using the calibration given by Tremonti et al. (2004). The small blue points refer to the observational data on metal-rich branch, and the 
stars refer to the SDSS sample galaxies with lower metallicitics, 12+log(O/H)=8.0-8.5, which arc the Bayesian metallicity estimates from 
Tremonti et al. (2004) and are in the "transition region" from the metal-poor to the metal-rich branches. The open triangles in cyan color 
represent the 67 H II regions in 21 dwarf irregular galaxies taken from van Zee & Haynes (2006). The "primary" (the dot-dashed line), 
"secondary" (the long-dashed line) components and the combination of these two components (the solid line) taken from Vila-Costas & 
Edmunds (1993) have also been plotted. The dotted line refers to the prediction of a numerical one-zone model from Henry et al. (2000). 
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Fig. 9. — log(N/0) abundance ratios of the SDSS sample galaxies as a function of the Bayesian metallicity estimates 12+log(O/H)T04 
from Tremonti et al. (2004). The lines are the same as in Fig. IS1 
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Fig. 10. — log(N/0) abundance residuals of our SDSS sample galaxies as a function of the equivalent widths of their Ha emission lines 
EWs(Ha) (in A). The solid line refers to a linear least-squares fit with a rms ~0.057 for the data points (Eq.(5)). 
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TABLE 1 

Coefficients for the oxygen abundance diagnostics. Defining oh=12+log(0/H) here; assuming y=b0+b1x, and y=a0+a1x+a2x 2 +a3x 3 for the 

LINEAR AND 3-ORDER POLYNOMIAL FITS, RESPECTIVELY. HERE, Y— LOGi?23 ; 12+L0G(0/H) R23 , 12+LOG(O/H) T0 4 and 12+LOG(0/H) K 99 FOR THE 

different cases, and x= N2, 03N2, N202, N2S2, 03Hb for. each of the "y" case, a is the rms deviation from the data to the fitted line in dex. 
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